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Abstract 

The interplay between impurities and interactions is studied in the gap- 
less phase of two-chain Hubbard model in order to see how the screening of 
impurity potentials due to repulsive interactions in single-chain model will 
be changed by increasing the number of channels. Renormalization group 
calculations show that charge stiffness, and hence persistent current, of the 
two-chain model are less enhanced by interactions than single chain case. 

I. INTRODUCTION 

Interacting electrons in a disordered system constitue one of the most challenging prob- 
lems. Questions on roles of interactions and disorder in quantum wires were again cast by the 
experiments of persistent current in mesoscopic rings (lj . The ground state of a ring system 
with a magnetic flux breaks time reversal symmerty and has a finite circulating current as a 
periodic function of the flux. This persistent current was predicted theoretically long time 
ago |2| , but we had to wait a few decades until technological progress allowed to observe the 
current. The experiments revealed that the ring supports a large current even with mod- 
est amount of impurities. Non-interacting electrons with impurities are shown || to have 
persistent current with a factor l/L multiplied, where I is the elastic mean free path and L 
the ring size. The current observed in the experiments was much larger than this value and 
other effects had to be taken into account to explain the large current. Electron-electron in- 
teractions were considered to be the first candidate as an origin of the large current, because 
a simple consideration gives a scenario that repulsive interactions would prevent particles 
from gathering in deep impurity holes and hence the interacting electrons would not be 
easily pinned by impurities, compared to non-interacting electrons. However the problem 
is not that simple as shown in the studies of one-dimensional (ID) system. For example, 
the roles of interactions are quite different between the systems with spinning and spinless 
electrons. The detail of the difference is summarized in section II. 

Since ID systems allow us to study the physical properties rather rigorously with the 
help of bosonization and renormalization group (RG) techniques. Approaching the problem 



of persistent current from a viewpoint of ID systems, hence, has this kind of advantage. 
Analytical and numerical studies have been done for interacting electrons in a single chain 
with disorder potentials. Summary of their results is given also in section II. 

In order to approach more realistic model, one need to consider multi-channel effects 
on the chain. Recent growing interests in few-coupled chains are mostly focused on pure 
systems and there have been only a few studies on ladder electron systems with impurity 
potentials. In this paper we study two-chain Hubbard model with impurities as a first step 
from a single chain towards a finite cross section ring. In section III we derive RG equations 
of the renormalized strengths of impurity potentials in the gapless phase, and show how the 
relation between interactions and disorder is changed by increasing the number of channels. 
Summary is given in section IV. 

II. SINGLE CHAIN 

Let us first begin by summarizing the impurity-interaction effects on persistent current 
of single chains. For spinless electron models with impurities, analytical and numerical 
studies showed that persistent current is suppressed by interparticle repulsive correlation as 
far as the impurity potentials are weak enough A RG calculation gives the effective 

backward impurity scattering W is renormalized as dW/dl = (3 — 2K p )W where K p is the 
Luttinger parameter. @] The larger repulsive interactions give rise to the smaller K p and 
therefore to the stronger effective impurity potentials, which would diminish the persistent 
current. This is due to the enhancement of charge density wave (CDW) correlations in the 
ground state of the system. The repulsive interactions would enhance the CDW correlations, 
which make the system easily pinned by the impurities therefore make the persistent current 
smaller than the value of non-interacting systems. 

In the presence of the strong impurity potentials, on the other hand, the wave functions 
are almost localized. When the interactions are turned on between the localized particles, 
they start to escape from each other and the localization of the particles in the deep impurity 
potentials is relaxed effectively. As a consequence, the interactions tend to enhance the 
persistent current of systems with strong impurity potentials ||. 

In the experimental point of view, the above results are disappointing because the ma- 
terials used in the experiments are all in the diffusive or ballistic regime so the scenario of 
strong impurity potentials cannot be applied. 

The situation changes drastically, however, when the spin degrees of freedom are taken 
into account |7|]. The renormalized impurity potential is given by dW/dl = (3 — K p — 
K„ — y)W, where K v {y = p,a) is the Luttinger parameters of charge and spin sectors 
respectively and y measures the backward scattering strength between electrons of opposite 
spins. For Hubbard model with small U, this RG equation becomes 

dW/dl = (1 - U/ttv f )W, (1) 

where Vp is the Fermi velocity. Here electron-electron interactions would screen impurity 
potentials. Since spin density wave (SDW) correlation is dominant in the ground state of 
the repulsive Hubbard model, the interaction U makes the particle density uniform, and 
therefore make the coupling of the density to the impurities weak. This role of interactions 
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is a good news to explain the experimental data: electron-electron interactions weaken 
the impurity effects and the persistent current would be enhanced from the value of non- 
interacting systems. 

Note that, while 2kp component of density- impurity coupling is renormalized to smaller 
value due to the interactions as shown, 4kp component grows in the presence of the inter- 
actions. Since kp of spinning electrons is equal to 2kp of spinless fermions, 4k p coupling in 
infinite U Hubbard model corresponds to 2kp coupling of spinless model where the 2k f cou- 
pling is enhanced by the interactions as stated above. Therefore as U increases the effective 
strength of impurity potentials decreases as long as U is small and then increases for large 
U because of the effect of 4k F component ||. 

Those stories are about short range interactions. Then, what about long-range Coulomb 
interactions? Since it is shown || that long-range interactions make 4fcp density correla- 
tions dominant, the interaction would enhance the effective strengh of impurity potentials 
and hence suppress persistent current. The suppression of persistent current by Coulomb 
interactions was observed in a numerical calculation . 



Those analyses have almost revealed the interplay between interactions and impurities 
in single chain systems. For more realistic discussions, however, one has to consider multi- 
channel effects because the real ring systems have finite cross sections. Although it is 
shown that the persistent current of single-chain Hubbard model would be enhanced by 
the presence of electron-electron interactions, it is not clear whether this effect still remains 
in multi-channel systems and, more precisely, whether this effect will become stronger or 
weaker as the number of channels increases. In order to see this, we study two-channel 
Hubbard model as a starting point toward multi-channel models. 

III. TWO-CHAIN LADDER 

In this section we study two-chain Hubbard model which is described by H = H + H int ; 

H = ~t J2 c Li £ Wl ~t±J2 C Ui c 2si + h.c, (2) 
a,s,i s,i 

Hint = U y]n a iin a ii, (3) 



(V) 



where a (= 1, 2) is chain index. H can be diagonalized by a 0)7r = (ci ±c 2 )/v^2, where o and 
7r denote the bonding and anti-bonding bands respectively. Using the standard bosonization 
scheme for fermions with spins, we can describe the system by four fields cj)*, <p°, <p°. Then, 
introducing a linear combination of the fields: <p v ± = (<f)° ± v2 where v = p,a, we have 
the Hamiltonian, H = Hi + H2 + H3, 

H x = £ / p[u ur K ur (nU ur ) 2 + ^(d x cf> vr ) 2 }, (4) 
v = p,a vr 
r — ± 



H 2 = 

2 na 



1 f 
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H 3 = J dx{gW 00 cos 20 ff _ + g^ Q cos 20 p _ + g^ Q7T cos 26 a _ - cos 26 p _} cos 20 CT+ , (6) 

where a is the lattice constant, if 2 is associated with interband processes induced by intra- 
chain forward scattering, and if 3 with the intrachain backward scattering. The sign of each 
term in Eqs.(|5|) and ((^) is determined according to a simple algebra of Majorana Fermions 
ljj , introduced to preserve the proper anticommutation relations between Fermion fields 
with differing band and spin indices. The interaction parameter g a ^s represents the scatter- 
ing from (5,7) to (a, (3). All g's in Eqs. (|5|) and (H) are equal to U for the case of Hubbard 
model. We used in Eqs. (^) and (|B|) that g^J 00 = g^} n7! . and the similar properties of g's. 
The impurity scatterings are described by 

if = V V imp n ■ (7) 

a,s,i 

Applying the bosonization to H imp and separate it into the forward and backward scattering 
parts, H imp = H f + H b , we get 

H f = J dxrn(x){d x (f> p + d x <fr p ) (8) 

+ / dx^[e i{<t '^ + ^- ) cos(0 CT _ + 6 a -) + e i(-*P-+»P-) cos(0 CT _ - 6 a _) + h.c], 
J na 

H b = — f dx£i(x){e- i « , ' ++ ^- ) cos(0 (J+ + <f> a J) + e -^+-^-) C os(0 (T+ - + h.c. 

na J 

+ — / dx&(x){e- i{ ^ ++(f> ° +) cos(V + d a ) + e -^ P +-^+) cos (6 - 6 a )} + h.c. (9) 
na J 

r)i and £1 are the random potential fields within o and n bands, and 772 and £2 are the random 
hopping between o and n bands. In order to make a quantitative comparison with the single 
chain case, we use the same distribution for rj and £ as used in the single chain model. In the 
single chain model we assumed Gaussian distribution and 77(0)77(2;) = D5(x), £*(0)£(x) = 
W5(x). We assume the same for V^ mp , that is, we have 774(0)77^(0;) = (Di/2)5ij5(x) and 
£*(0)£ 3 -(ar) = (Wi/2)8ijS(x), where i, j = 1, 2. The factor 1/2 appeared in the transformation 
from ci >2 to a 0)7r . 

The Hamiltonian H has been studied analytically []XTj — P~^] and numerically [T5-(H| and is 
shown to have a rich phase diagram. The model including impurities, that is H + H imp , was 
recently investigated using a RG method pOyZIfl . It was found in a spin gap phase of the 
repulsive Hubbard ladder that the dominant coupling between charge density and impurity 
potentials is the 4kp Fourier component, which is already shown in the single chain models 
to be reduced by the renormalizations. Therefore the persistent current of the spin gap 
phase is suppressed by the repulsive interactions, at variance with the single Hubbard chain. 

The spin gap phase has attracted a lot of attentions in connection with the possibility 
of superconductivity. But here our final concern is to see the interplay between interactions 
and impurities in multi-channel systems, keeping the experimentally relevant situations in 
mind. Since the ground states are gapless in the two limits, a single channel and the real 
ring with many channels, we would like to focus on the gapless phase of the two-chain model 
|22[ in order to start a systematic interpolation between the limits, although the gapless 



phase occupies only a small area in the phase diagram. [13 
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Here we first study the effect of term. Using the replica trick for the random variables 
and following the standard RG procedure, we get the RG equations of the backward impurity 
scatterings W\ and W 2 in the gapless phase; 



dWi 
dW 2 



[3 - (K p+ + K p „ + K a+ + K 9 J)/2 - (jW + g%J/2™*]Wi 
[3 - (K p+ + ^L + K(T+ + 1 r^)/2- (g^ - g^)/2-Ku a }W 2 . 



(10) 
(11) 



The renormalizations to g^> in Eqs. ( [XQ|) and (|11]) are described by 
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dl 


= (2- 
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(12) 
(13) 
(14) 
(15) 



In the case of Hubbard model with small U and small disorder, we ignore the renormal- 
ization to K and g and then we get 



dWi 

~df 

dW 2 

~df 



[1 - U/nv F ]W!, 
[l + 0((U/v F ) 2 )]W 2 



(16) 
(17) 



The coefficient of (U/v F ) 2 term in Eq. ( |TT| ) is positive. First, we notice that Eq. ([R]) has the 
same form as Eq. (|]), and the renormalization to the backward impurity scatterings within 
the bands is unchanged even when two chains are coupled, which is to be contrasted to the 
case of spin gap phase P0 |. The impurity scatterings between the bands behave differently, 
however. The presence of the electron-electron interactions leads to the larger W 2 , although 
the dependence of W 2 on U is rather weak since the leading order in Eq. ( |iT|) is (U/v F ) 2 . 

The current coupled to the flux penetrating the ring, is J2s(hs + 32s) = J2 s (jos + j-ns) = 
(2/ y / 7r)n p+ . Therefore only (p+) mode contributes to persistent current and charge stiffness 
D, where the latter is given by 4K p+ u p+ . Then the RG equation of charge stiffness is 



dD 
~df 



2DK p+ u p+ a 



(18) 



where W + = (Wi + W 2 )/2. Remember W\ behaves the same way as in the single-chain case, 
namely U suppresses W\, whereas W 2 increases with U. The RG equation of D for single 
Hubbard chain is dD /dl = — (2DK p u p a/7ru^)W. One then see D of the gapless phase of two 
chain model is enhanced by the interaction U, which is in contrast to the spin gap phase, but 
the enhancement is weaker than single chain. Namely, the enhancement of charge stiffness 
(therefore persistent current) due to the electron-electron repulsive interactions becomes 
smaller by doubling the number of channels. 
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The forward scatterings from impurities, Hf, were already studied in Ref. [fj(],^l[]. The 



first term in Hf can be absorbed in transforming the definition of (p p ± and then has no 
effect. The second term, however, was shown to make the interaction terms effectively 
weak, which results in that the enhancement of the persistent current due to the presence 
of the interactions would become even smaller in two chains than in single chain. 



IV. SUMMARY 



We studied the impurity effect on the gapless phase of the two-chain Hubbard model, as 
a first step from a single chain analysis toward the finite cross section ring. In contrast to 
the spin gap phase, the RG calculations show that the charge stiffness, and hence persistent 
current, would be enhanced by the repulsive interactions. The enhancement, however, is 
smaller than in single Hubbard chain. This indicates a possibility that the large persistent 
current observed in the real rings with finite cross sections might not be explained by simply 
extending the mechanism of interaction-enhanced persistent current in single Hubbard chain. 
In order to draw more definite conclusion, one has to check on three chain system and work 
in this direction is in progress. 
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